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ABSTRACT 

The jet opening angle and the bulk Lorentz factor Fq are crucial parameters for the com- 
putation of the energetics of Gamma Ray Bursts (GRBs). From the ^30 GRBs with measured 
6'jct or Fq it is known that: (i) the real energetic E^, obtained by correcting the isotropic equiv- 
alent energy Eiso for the collimation factor ^ 6*?,^, is clustered around 10^°-10^^ erg and it 
is correlated with the peak energy E-^ of the prompt emission and (ii) the comoving frame 
i?p and E'^ are clustered around typical values. Current estimates of Fq and 6'jct are based on 
incomplete data samples and their observed distributions could be subject to biases. Through 
a population synthesis code we investigate whether different assumed intrinsic distributions 
of Fo and 6'jct can reproduce a set of observational constraints Assuming that all bursts have 
the same E'^ and E'^ in the comoving frame, we find that Fq and cannot be distributed 
as single power-laws. The best agreement between our simulation and the available data is 
obtained assuming (a) log-normal distributions for 0jet and Fq and (b) an intrinsic relation 
between the peak values of their distributions, i.e 6'jct^'^Fo=const. On average, larger values 
of Fq (i.e. the "faster" bursts) correspond to smaller values of 6'jct (i-C- the "narrower"). We 
predict that ^6% of the bursts that point to us should not show any jet break in their afterglow 
light curve since they have sin 6'jct < l/ro- Finally, we estimate that the local rate of GRBs 
is ^0.3% of all local SNlb/c and ^4.3% of local hypernovae, i.e. SNIb/c with broad-lines. 
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1 INTRODUCTION 

Gamma Ray Bursts (GRBs) have extremely high energetics. The 
isotropic equivalent energy -Biso, released during the prompt phase, 
is distributed over four orders of magnitudes in the range lo^''"^* 
erg. Eiso correlates with E^, i.e. the peak of the vFi, spectrum 
(Amati et al. 2002, 2009): Ep oc . This holds for long duration 
GRBs. A similar correlation exists between the isotropic equiva- 
lent luminosity Liso and Ep (Yonetoku et al. 2004) obeyed also by 
short events (Ghirlanda et al. 2009). The scatter of the data points 
around the _Ep — Eiso correlation, modeled with a Gaussian, has a 
dispersion Gbc ~ 0.23 dex (see e.g. Nava et al. 2012 for a recent 
update of these correlations). This dispersion is much larger than 
the average statistical error (t_Bj^^ = 0.06 dex and (t_e ^j,,^ = 0.10 
dex associated with _Biso and Ep, respectively. 

Since _Biso is computed assuming that GRBs emit isotropi- 
cally, it is only a proxy of the real GRB energetic. GRBs are thought 
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to emit their radiation within a jet of opening angle 6jct ■ If the jet 
opening angle 9jet is known, the true energy E^c^EisoOf^t and the 
true GRB rate can be estimated (Frail et al. 2001). 

The estimate of 9 jet is made possible by the measure of the jet 
break time f break, typically observed between 0.1 to >10 days in 
the afterglow optical light curve. Although Sjct has been measured 
only for ~30 GRBs (Ghirlanda et al. 2007) it shows that: 

(i) E^ clusters around 10^" erg with a small dispersion (Frail et 
al. 2001; but see Racusin et al. 2009; Kocevski & Buder 2008); 

(ii) E^ is tightly correlated with Ep (Ghirlanda, Ghisellini & 
Lazzati 2004; Ghirlanda et al. 2007) with a scatter CTsc = 0.07 dex 
(consistent with the average statistical error aE^, = <^E-p — 0.1 dex 
associated with E^, and Ep)\ 

(iii) the true rate of local GRBs ranges from ~ 250 Gpc^'^ yr~^ 
(e.g. Frail et al. 2001) to ~ 33 Gpc"^ yr"^ (Guetta, Piran & Wax- 
man 2005). These different values are mainly due to the different 
values assumed for the collimation factor / oc 97^ . The true GRB 
rate can be compared with the local rate of SN Ib/c (e.g. Soderberg 
2006; Guetta & Delia Valle 2007; Grieco et al. 2012), i.e. the can- 
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didate progenitors of long GRBs, and allows to estimate the rate of 
orphan afterglows (e.g. Guetta et al. 2005). 

The Ep — Eiso, Ep — Lisa and Ep — Ej correlations could 
enclose some underlying feature of the ORB emission mecha- 
nism (e.g. Rees & Meszaros 2005; Ryde et al. 2006; Thompson 
2006; Giannios & Spruit 2007; Thompson, Meszaros & Rees 2007; 
Panaitescu 2009), of the GRB jet structure (e.g. Yamazaki, loka 
& Nakamura 2004; Eichler & Levinson 2005; Lamb, Donaghy 
& Graziani 2005; Levinson & Eichler 2005) or of the progenitor 
(e.g. Lazzati, Morsony & Begelman 201 1). An intriguing appUca- 
tion of these correlations is the use of GRBs as standard candles 
(Ghirlanda, Ghisellini & Firmani 2005; Firmani et al. 2005; Amati 
et al. 2009). 

The presence of outliers of the Ep — Eiso correlation in the 
CG/?0/BATSE GRB population (Band & Preece 2005; Nakar & 
Piran 2005; Shahmoradi & Nemiroff 2011) and in the FermilGBM 

burst sample (Collazzi et al. 2012) and the presence of possible in- 
strumental biases (Butler et al. 2007; Butler, Kocevski & Bloom 
2009; Kocevski 2012) caution about the use of these correlations 
either for deepening into the physics of GRBs and for cosmologi- 
cal purposes. Although instrumental selection effects are present, it 
seems that they cannot produce the correlations we see (Ghirlanda 
et al. 2008; Nava et al. 2008; Ghklanda et al. 2012b). Moreover, a 
correlation between Ep and Liso is present within individual GRBs 
as a function of time (Firmani et al. 2009; Ghirlanda et al. 2010; 
2011; 201 la), suggesting that the radiative process(es) might be the 
origin of the Ep — Lisa correlation. Despite these studies, the spec- 
tral energy correlations of GRBs and their possible applications are 
still a matter of intense debate. 

A new piece of information recently added to the puzzle is that 
the GRB energetics (E{g„, L{g^ and Ep) appear nearly similar in the 
comoving frame (Ghirlanda et al. 2012 - G12 hereafter). To mea- 
sure these comoving quantities^ we have to know the bulk Lorentz 
factor To, that can be estimated through the measurement of the 
peak time tpeak of the afterglow light curve. G12 could estimate 
To in 30 long GRBs with known z and well defined energetics, 
finding that: 

(i) £;iso(iiso)oc and BpOcLo; 

(ii) the comoving frame £i'so~3.5x 10®^ erg (dispersion 0.45 
dex), L[^r^5xW^^ erg s"^ (dispersion 0.23 dex) and Sp~6 keV 
(dispersion 0.27 dex). 

These results imply that the Ep — Eiso and Ep — Liso correlation 
are a sequence of different To factors (see also Dado, Dar & De 
Rujula 2007). 

The Sjet values of GRBs are known only for a couple of 
dozens of bursts (Ghirlanda et al. 2007). 6'jet appears distributed 
as a log-normal with a typical ^jot~ 3° (Ghirlanda et al. 2005). 
By correcting the isotropic comoving frame energy E'i^o by this 
typical jet opening angle, the comoving frame true energy E!y re- 
sults ~ 5 X 10** erg. In G12 we also argued that in order to have 
consistency between the Ep — Ej and the Ep — Eiso correlations 
one must require ^j^etTo = constant. A possible anti-correlation be- 
tween 6 jet and Fq is predicted by models of magnetically acceler- 
ated jets (Tchekhovskoy, McKiimey & Narayan 2009; Komissarov, 
VlahaMs & Koenigl 2010) but, at present, only 4 GRBs have an 
estimate of 6'jot and Fq and well constrained spectral properties. 

The measure of ^jct relies on the measure of tbroak, that in 
turn requires the follow up of the optical afterglow emission up to 

1 Primed quantities are in the comoving frame of the source. 



a few days after the burst explosion (Ghirlanda et al. 2007). The 
measurement of tbreak is difficult, not only because it requires a 
large investment of telescope time, but also because several tbrcak 
are chromatic (contrary to what predicted; but see Gliisellini et al. 
2009), and the jet break can be a smooth transition whose measure- 
ment requires an excellent sampling of the afterglow light curve 
(e.g. Van Eerten et al. 2010, 2011). Another complication is that 
the early afterglow emission is characterized by several breaks. 
For instance, the end of the plateaux phase typically observed in 
the X-ray light curves, if misinterpreted as a jet break, biases the 
^jct distribution towards small values of Ojct (Nava, Ghisellini & 
Ghirlanda 2006). Finally, the measure of large 9 jet is complicated 
by the faintness of the afterglow and its possible contamination by 
the host galaxy emission and the supernova associated to the burst. 
Several observational biases could shape the observed ^jct distri- 
bution. Among these the fact that more luminous bursts (i.e. those 
more easily detected) should have the smallest jet opening angles. 
For all these reasons the observed distribution of ^jet might not be 
representative of the real distribution of GRBs jet opening angles. 

The distribution of Fq is centered around Fo=65 (130) in 
the case of a wind (uniform) density distribution of the circum- 
burst medium. The distribution of Fq is broad and extends between 
Fo~20 and Fo^SOO. These results are still based on a sample of 
only 30 GRBs (G12). The difficulties of early follow-up of the op- 
tical afterglow emission could prevent the measure of very large Fq 
on the one hand, while the possible contamination by flares (Bur- 
rows et al. 2005; Falcone et al. 2007) or by other (non afterglow) 
emission components (e.g. Ghisellini et al. 2010) at intermediate 
times could prevent the estimate of the low-end of the Fo distribu- 
tion. One could argue if GRBs can have Fq of a few. While there 
are some hints that GRB060218 should have Fo~5 (Ghisellini et al. 
2006) the classical compactness argument, for typical GRB param- 
eters (e.g. Piran 1999), requires that Fo^ 100-200. This argument 
was successfully applied to few bursts observed up to GeV ener- 
gies by LAT on board Fermi (e.g. Abdo et al. 2009, Ghirlanda et al. 
2009) to derive lower limits of several hundreds on Fq. If, instead, 
the highest energy photon detected has an energy of say -Emax ~ 3 
MeV, the lower limit derived from the classical compactness argu- 
ment would be ro> a few (i.e. ~ 2Eui^^/rne(?). Therefore, also 
in the case of Fo, the observed distribution, derived with still few 
events, could be not representative of the real distribution of this 
parameter. 

The main aim of this paper is to constrain the distribution of 
Fo and ^jet in GRBs using the available independent constraints. 
This aim can be translated into a simple question: do Ojct and Fo 
follow power law distributions or do they follow some kind of 
peaked distribution (e.g. a broken power law or a log-normal)? In 
both cases the resulting distributions could be different from the ob- 
served ones since some selection effect (as discussed above) might 
prevent to measure very low and/or high values of Fo and 6'jct. An- 
other scope of the present paper is to test which is (if any) the rela- 
tion between Ojet and Fo. A relation ^jet^Fo=const was assumed in 
G12 to explain the spectral energy correlations and a similar rela- 
tion seems to arise from numerical simulations of jet accelerations 
(Tcheckolskoy et al. 2012). Here we use several observational con- 
straints and test whether there is a ^jct'^Fo=const relation and try 
to constrain its exponent a. One important effect that we consider 
in this paper for the first time is the colUmation of the burst radi- 
ation when Fo is small. In general we are led to think that given 
a value of the collimation corrected energy E~,, the correspond- 
ing isotropic equivalent energy is Eiso^Ej/9jet^ . This is true if 
the beaming of the radiation is "dominated" by the jet opening an- 



© 0000 RAS, MNRAS 000, 000-000 



Characteristic 6'jet and Tq in GRBs 3 



gle, i.e. l/Fo^J sin6'jct. However, GRBs with very low To could 
have l/ro> sin^jct and in this case the isotropic eqtiivalent en- 
ergy is determined by To (i.e. Eiso^E-iVo^ - see §. 2.2) rather than 
by ^jet. This effect, introduces a limit (oc^iso^''^) in the classical 
i?p — E\so plane (Fig.l) accounting for the absence of bursts with 
intermediate/low Bp and large values of i?i,,o. This limit can also 
partly account for the problem of "missing jet breaks" since these 
bursts with l/Fo ^ sin^jet should not show any jet break in their 
afterglow light curve (§4.5). 

We rely on a GRB population synthesis code that we have re- 
cently adopted to explore the issue of instrumental selection biases 
on the Ep — Liso correlation (Ghirlanda et al. 2012b). 

The simulation steps are described in §2 while the observa- 
tional constraints that we aim to reproduce are outlined in §3. In 
§4 we present our results. We summarize and discuss our find- 
ings in §5. Throughout the paper a standard flat universe with 
/i = = 0.7 is assumed. 



2 POPULATION SYNTHESIS CODE 

So far, the approach adopted in studying the spectral-energy corre- 
lations and the distributions of ^jct or To was (i) to derive the col- 
limation corrected E-p — E^ correlation by correcting the isotropic 
energy Eiso for the coUimation factor oc ^?et (e-g- Ghirlanda et al. 
2004), or (ii) to derive the comoving frame properties of GRBs by 
correcting, for the Fq factor, the isotropic values Eiso, Liso and Ep 
(G12). 

In this paper, we tackle the problem from the opposite side and 
jointly work with ^jet and Fq: we assume that GRBs have all the 
same comoving frame E'p and E!y and simulate GRB samples with 
different distributions of Fo and 6'jct. This produces a population of 
GRBs with known energetics iJiso, peak energy Ep and observer 
frame fluence F and peak flux P. We would like to stress that our 
main assumption (same E'p and E!y for all burst) is a crude simpU- 
flcation. Nevertheless, our assumption can work if the real Ep and 
E!y distributions are indeed narrower that the distributions of the 
corresponding observed quantities. Recently, Giannios 2012 have 
shown that in photospheric models a comoving frame peak energy 
E'pr^l.5 keV is expected. 

The observational constraints that we aim to reproduce (see 
§3) are: (i) the rate of GRBs observed by Swift/BAT, CGRO/BATSE 
and FermilGBM, (ii) the Ep — Eiso correlation defined by the com- 
plete sample of Swift bright bursts (Salvaterra et al. 2012; Nava et 
al. 2012) and (iii) the fluence and peak flux distributions of the pop- 
ulation of bursts detected by Fermi/GBM (Goldstein et al. 2012) 
and CGRO/BATSE (Meegan et al. 1998). 

Note that, since one of the aims of the present paper is to con- 
strain the distributions of Sjct and Fo we cannot adopt the observed 
ones (discussed in the introduction) as constraints, otherwise we 
would fall into a circular argument. The distributions of Fo and 9jet 
that we assume in our simulations (power law, broken power law, 
log-normal) have all their characteristic parameters (slope, normal- 
ization, break values, width etc.) free to vary. These parameters are 
what we aim to constrain through our population synthesis code. 

In Fig. 1 we show the rest frame peak energy Ep versus the 
total energy E (where E here is generically used to indicate an 
energy, either isotropic or collimation corrected). We highlight dif- 
ferent regions (I, II and III) that arc useful to explain the simulation 
steps (§2. 1). This plane will be one of our observational constraints: 
in Fig. 1 we show (black filled points) the Swift complete sample 



of bursts (Salvaterra et al. 2012; Nava et al. 2012) which we aim to 
reproduce through our simulations. 

2.1 Simulation steps 

Our starting assumption is that all GRBs have the same comov- 
ing frame E'p=l.5 keV and E!y=1.5xlQ'^^ erg. This is shown by 
the black circle in Fig. 1. G12 find that Lj^o^const and that the 
observed duration Too does not depend on Fo. Therefore, in the co- 
moving frame, Tgo oc FoTgo oc Fq. It follows that E'^=L'isoTgo9fet 
is also constant if, as discussed in G12, 9?3tFo=const. Although 
some dispersion of the values of Ep is present in the sample of 
G12, the value of E'p that we assume here is consistent at the la 
level of confidence with the distribution of values reported in G12 
for the wind density ISM. 

The main steps of our simulation are: 

(i) we simulate a population of GRBs distributed in redshift z 
between z = and z = 10 according to the GRB formation rate 
(GRBFR) ip{z). This is formed by two parts: ip{z) = e{z)R{z). 
The first term is a cosmic evolution term, while R{z) is taken from 
Li (2008) (which extended to higher redshifts the results of Hopkins 
& Beacom 2008): 

_ 0.0157 + 0.1180 
' 1 + («/3.23)4-66 

R{z) is in units of Mq yr~^ Mpc~^. Concerning e{z), Salvaterra 
et al. (2012) derived the luminosity function of GRBs by jointly fit- 
ting the redshift distribution of a complete sample of bright GRBs 
detected by Swift and the count distribution of a larger sample of 
BATSE bursts. They found that either the evolution of the lumi- 
nosity function or the evolution of the density of GRBs is required 
in order to account for these data sets. We assimie the same term 
e{z) = (l + a)^ '' found by S12. 

(ii) We assign to each GRB a bulk Lorentz factor Fq extracted 
from a specified distribution, in the range [1, 8000]. The upper limit 
(Fo,max = 8000) is somewhat arbitrary, but large enough to en- 
compass all the values of Fq estimated so far, and in particular the 
large values derived for the few GRBs detected by the LAT instru- 
ment on board Fermi, if the GeV emission is interpreted as after- 
glow (Ghisellini et al. 2010). 

For each simulated burst the rest frame peak energy Ep and the 
energy E-^ are (see G12): 

5Fo 



-^peak — -^peak f 



where Fo=l/(l — /3o)^^^ - The simulated bursts define a correlation 



2/3o 



E'^Fo 



(2) 



between Ep and Ey. 



E'^ 



peak 



peak 



E' 



(3) 



5-2A) 5-2/3o 
for /3o ~ 1 this corresponds to the Ep — Ej correlation in the 
case of a wind density profile (Nava et al. 2006). This relation is 
shown in Fig. 1 with the solid black line (labelled EpOcE^). The 
simulated Fq distribute Ep between 1.5 keV (Fo=l) and ~20 MeV 

(r()=8ooo). 

(iii) We assign to each simulated burst a jet opening angle 6'jotG 
[1°, 90°] extracted from a specified distribution. 

(iv) The probability for a burst to be observed from the Earth 
depends on the viewing angle 6'view between the jet axis and the 
line of sight of the observer. We extract randomly a viewing angle 
Sview from the cumulative distribution of the probability density 
function sin^view 
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Figure 1. Rest frame plane of GRB energetics. The large black dot corresponds to the main assumption of our simulations, i.e. that all bursts have similar 
comoving frame i?p=1.5 keV and i?^^=1.5 X 10'** erg. Assigning a certain Fq to the burst, this moves along the line E-p oc Ej. Since Fq > 1 and we assume 
a maximum Fq of 8000, regions (I) are forbidden. Since all our simulated bursts have 9jt,t^90°, they cannot lie in region (II). When Fq is small, the beaming 
cone ~ 1 /Fq can become wider than the aperture of the jet. In this case the isotropic equivalent energy becomes E{so=E^ (1 + Po)^'q, that is smaller than 

1 /3 

the energy calculated through Eiso=i?-y/(l-cos9jet). This introduces a limit E-p oc E.^^^ and bursts cannot lie on the right of this limit. Consequently, region 
(III) is forbidden. The black dots con'epond to the real GRBs of the Swift complete sample. The fit to the Swift complete sample is shown by the dot-dashed 
line. 



(v) In order to compare the simulated bursts with the source 
count distribution of existing samples of GRBs (see §3) we com- 
pute the observer frame peak fluxes P and fluences F. To this 
aim we assume a typical spectrum described by the Band func- 
tion (Band et al. 1993), with low and high photon spectral indexes 
Q = —1.0 and /3 = —2.3, respectively (i.e. corresponding to the 
typical values observed by different instruments - e.g. Kaneko et 
al. 2006; Sakamoto et al. 2011)^. The fluence F of each simulated 
burst in a given energy range is computed by re-normalizing this 
spectrum through the bolometric fluence -Fboi=£'iso(l + z)/47rdL, 
where is the luminosity distance for a given redshift z. To de- 
rive the peak flux P, we assign to each burst an (observer frame) 
duration Tgo extracted from a distribution centered at 27.5 s and 
with a dispersion crLogTgo = 0.35. This distribution is truncated 
at Tgo = 2 s because we consider only long duration GRBs in 
this analysis. Such a duration distribution is similar to that of the 
FermilGBM GRBs (Paciesas et al. 2012; Goldstein et al. 2012) and 
includes also very long bursts with Tgo ~300 s. We assume that the 
bursts have a simple triangular light curve and derive the peak lu- 
minosity as Lpoak = 2i?iso(l + z) /Tgo. The peak flux P in a given 
energy range is obtained by re-normalizing the spectrum through 
the bolometric peak flux Pboi = I/pea,k/47rdL- 



^ These values are also assumed by S 12 to constrain the LP of GRBs. 



2.2 Computation of Eiso 

The isotropic equivalent energy Eiso of the simulated bursts can 
be derived from E^. Since Sjot^90°, simulated bursts cannot be in 
region II of Fig. 1 and i^iso can take values on the right hand side 
of the limit of Eq. 3 shown in Fig. 1 . According to the values of 9 jet 
and To assigned to each simulated bursts, the isotropic equivalent 
energy is: 

Siso = E^/{1 - cosejet) if l/Fo Sin6ljet (4) 

Bi.o = E^{1 + if 1/ro > sinSjct (5) 

In the latter case E^so is smaller than in Eq. 4. This introduces a 
limit in the Ep — Eiso plane of Fig. 1 corresponding to the line: 

E- 1 ^''^ 

^^^^^ i(5-2/3o)^(l + /?o)J 

(labelled £'pOc£iso^''^ in Fig. 1). For a given 6'jct, bursts with a 
small value of To will have an Eiso computed through Eq. 5 and 
will lie on the limiting line of region III in Fig. 1 . Their radiation 
is, indeed, collimated within an angle arcsin(l/ro) which is larger 
than their 6'jot. 

Simulated bursts can populate the region delimited by bound- 
aries (I, II and III) in Fig. 1. This is one (among others) observa- 
tional constraint that we will adopt in our simulations (§3) to con- 
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strain the distributions of To and ^jct and tiieir possible relation. 
According to the relative values of 6 jet, To and 6view, simulated 
bursts are classified as: bursts "pointing to us" (PO, hereafter), i.e. 
those that can be seen from the Earth, with sin^?vicw^max[sin^jet, 
l/Fo] and bursts pointing in other directions (NPO, hereafter), i.e. 
not observable from the Earth, with sinSviow>max[sin6jot, l/Fo]. 
We will compare the PO simulated bursts with our observational 
constraints, while the entire population of simulated bursts (i.e. PO 
and NPO) will be used to infer the properties of GRBs (e.g. the 
distributions of To and Sjet and the true burst rate). 



3 OBSERVATIONAL CONSTRAINTS 

In order to test whether ^jct and To assume characteristic values or 
not we compare the population of simulated bursts with real sam- 
ples of GRBs. In this section we describe our observational con- 
straints. We consider the ensemble of GRBs detected by the Burst 
Alert Telescope (BAT) on board Swift, the Gamma Burst Monitor 
(GBM) on board Fermi and the Burst And Transient Source Exper- 
iment (BATSE) on board the Compton Gamma Ray Observatory 
(CGRO). 

3.1 The Swift BAT complete sample 

Salvaterra et al. (2012 - S12 hereafter) constructed a sample of 
bright Swift bursts consisting of 58 GRBs detected by SwiftfBAT 
with P ^ Piim=2.6 ph cm"^ s"^ (integrated in the 15-150 keV 
energy range). Fifty four of these events have a measured redshift z 
so that the S 12 sample is 90% complete in redshift. Forty six (out of 
54) GRBs in this sample have well determined spectral properties 
(filled circles in Fig. 1) and define a statistically robust Ep — Eiso 
correlation with rank correlation coefficient p = 0.76 and chance 
probability P = 7 x 10"^° (Nava et al. 2012, N12 hereafter)^. 
The correlation properties (slope and normalization) of the com- 
plete Swift sample are consistent with those defined with the incom- 
plete larger sample of 136 bursts with known z and spectral param- 
eters (see N12). Therefore, the distribution of the Swift complete 
sample (46/54 events with well constrained Ep) in the Ep — Eiso 
plane is representative of the larger (heterogeneous) population of 
GRBs with measured z and well constrained spectral properties. 
The 46 GRBs of the complete Swift sample define a correlation 
Ep oc ^w-eiio oe (shown by the dot-dashed line in Fig. 1) with a 
scatter (computed perpendicular to the best fit line) with a Gaussian 
dispersion a =0.29 dex. 

The Swift complete sample of S12 contains ~l/3 of the bursts 
detected by Swift^ with P ^2.6 ph cm"^ s"\ We verified that the 
Swift complete sample of 54 events selected by S 12 is representa- 
tive of the larger population of 149 long Swift bursts with P ^Pum- 
the Kolmogorov-Smirnov test on the peak flux distribution of the 
two samples gives a probability of 0.6 that the two distributions are 
drawn from the same parent population. These bursts were not in- 
cluded in the selection of S12 because they do not have favorable 
conditions for ground-based follow up. 

These 149 events with P ^Pum are the bursts detected by 
Swift in ~7 yrs from its launch within the (half coded) field of view 

^ The 8 GRBs without a secure estimate of the redshift or with incomplete 
spectral informations are consistent with the Ep — Eiso correlation defined 
by the 46 GRBs discussed here, see N12 for details. 
http://swift.gsfc.nasa.gov/docs/swift/archive/grb.table/ 



of ^ 1.4 sr of BAT. This corresponds to an average Swift detection 
rate of 7?.swift ~15 events sr^^ yr^- 

3.2 The Fermi GBM sample 

Another observational constraint that we consider is the population 
of bursts detected by the GBM on board Fermi. The spectral prop- 
erties of GBM bursts have been studied in Nava et al. (2011a) and 
compared to those of BATSE bursts in Nava et al. (2011b). More 
recently, the first release of the GBM spectral catalog (Goldstein 
et al. 2012) provided the spectral parameters and derived quanti- 
ties (i.e. peak fluxes and fluences) for 487 GRBs detected by the 
GBM in its first 2 years of activity. 398 bursts in this catalog are 
long events and have measured peak flux P and fluence F (both 
integrated in the 10 keV-1 MeV energy range)^. 

We cut the GBM sample to P > Pum = 2.5 ph cm"^ s"\ in 
order to account for the possible incompleteness of the sample at 
lower fluxes, obtaining 312 GBM bursts. 

The GBM is an all sky monitor that observes on average 
~60-70% of the sky. Therefore, the average GBM detection rate 
is 7?.GBM ~21 events sr^^ yr~^ with peak flux, integrated in the 
10 keV-1 MeV energy range, P ^ 2.5 ph cm~^ s~^. 



3.3 The CGRO BATSE sample 

We also consider the sample of GRBs detected by BATSE. The 4B 
sample (Meegan et al. 1998) contains 1540 long events and 1496 
of these have their P and F (both integrated in the 50-300 keV 
energy range) measured. The sample of 1496 BATSE bursts is cut 
at P ^ Piim = 1 ph cm"^ s"^ with 716 BATSE bursts above 
this threshold. Considering the average portion of the sky observed 
by BATSE, i.e. ~70% of the sky, the detection rate of BATSE is 
7?-BATSE ~16 events sr^^ yr~^ for GRBs with a peak flux, inte- 
grated in the 50-300 keV energy range, P ^ 1 ph cm~^ s~^. 

The lower detection rate of BATSE with respect to GBM is 
due to the different energy range where the peak fluxes are calcu- 
lated (i.e. 10 keV-1 MeV for GBM and 50-300 keV for BATSE, 
respectively). We verified that by considering the GBM bursts with 
peak flux P integrated in the same energy range of BATSE (i.e. 
50-300 keV) larger than 1 ph cm~^ s~^ (i.e. the same threshold 
adopted for BATSE), the GBM rate is equal to the BATSE one. 



3.4 Extraction of results 

From each simulation we extract three populations of GRBs among 
the bursts pointing to us (PO): 

(i) the Swift comparison sample: simulated GRBs with peak 
flux, integrated in the 15-150 keV band, larger than 2.6 ph cm~^ 
s^ ^ . We also require that their observer frame peak energy is in the 
range 15 keV-2 MeV. Indeed, this is the energy range where E"^" 
can be measured by presently flying satellites like Swift, Konus and 
Fermi. 

(ii) the GBM comparison sample: simulated bursts with a peak 
flux, integrated in the 10 keV-1 MeV energy range, larger than 2.5 
ph cm~^ s~^; 



P and F are reported in Goldstein et al. (2012) and were obtained by 
integrating the model that best fits the peak time resolved spectrum and the 
time averaged spectrum, respectively. 
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(iii) the BATSE comparison sample: simulated bursts with a 
pealc flux, integrated in the 50-300 IceV energy range, larger than 
l.Ophcm"^ s"\ 

The simulation is adjusted so that the Swift comparison sam- 
ple contains 149 GRBs, i.e. the same number of bright bursts de- 
tected by Swift (§3.1). Therefore, the Swift rate TZs^rUt is imposed. 
What we derive instead from the simulation is the rate of GBM 
and BATSE GRBs that we compare with the real rates of these two 
instruments described in §3.2 and §3.3 respectively. 

We also require that the Swift comparison sample is consistent 
with the Swift complete sample of S12. To this aim we compare 
them in the rest frame Ep — Eiso plane and in the observer frame 
^peJk ~ P plane deriving a 2 dimensional Kolmogorov-Smimov 
(KS) probability (one for the Ep - Ei^o and one for - F 

plane). Wc also verify through a 1 dimensional KS test that the red- 
shift distribution of the Swift comparison sample is consistent with 
that of the Swift complete sample. Finally we compare, through a 
ID-KS test, the fluence and peak flux distributions of the GBM 
and BATSE comparison samples with those of the real samples of 
GRBs detected by these instruments and described in §3.2 and §3.3, 
respectively. 

Since the Swift complete sample contains only the brightest 
Swift bursts it maps the high P end of the peak flux distribution of 
GRBs. The GBM and BATSE samples that we adopt here extend 
the comparison sample to lower values of P and ensures that our 
simulations reproduce also the faint end of the GRB population*^ . 

For each simulation we derive the following probabilities: 

• the 2D-KS probability that the Swift comparison sample is 
consistent with the complete Swift sample of S12 in the Ep — Eiso 
plane; 

• the 2D-KS probability that the Swift comparison sample is 
consistent with the complete Swift sample of S12 in the Ep^lj^ — F 
plane; 

• the ID-KS probability that the Swift comparison sample has a 
redshift distribution consistent with that of the S12 Swift sample; 

• the ID-KS probabilities that the GBM comparison sample is 
consistent with the GBM sample in terms of peak flux P and flu- 
ence F; 

• the ID-KS probabilities that the BATSE comparison sample 
is consistent with the BATSE sample in terms of peak flux P and 
fluence F; 

• we verify if the GBM rate predicted by the simulation is con- 
sistent, at la, with the GBM rate 7?.gbm. 

• we verify if the BATSE rate predicted by the simulation is 

consistent, at la, with the BATSE rate 7?-batse. 

For the KS probabilities we set a limit of 10~^ below which we 
consider that two distributions (either ID or 2D) are inconsistent at 
more than 3a. Each simulation, with its assumptions on the distri- 
bution of Sjet and To, is repeated 1000 times and we compute the 
percentage V of repeated simulations that produce GRB samples 
(i.e. Swift, GBM and BATSE comparison samples) consistent with 
our observational constraints. 



° The P values of the GBM sample are computed on the broad 10 keV-1 
MeV energy range (i.e. much broader than the 15-150 keV energy range 
of Swift). This ensures that the selected sample of the GBM bursts extends 
the population of GRBs to lower fluxes than those of the Swift complete 
sample. 



4 RESULTS 

In the following sections we present the results obtained with dif- 
ferent possible assumptions for the distributions of Fo and 6 jet- We 
want to test which one among the possible intrinsic distributions of 
To and Sjet that one can think of (e.g. power laws, broken power 
laws or log-normal) best reproduces the observational constraints 
described in the previous section. 

4.1 Power law distributions of Fo and Ojet 

We assume that both Ojet and Fo are distributed as power laws: 
dN/d9 jet oc djet and dN/dVo oc Fg. This corresponds to the hy- 
pothesis that Ojet and Fo do not have a characteristic value. We 
consider a £ [—2, —1] and c £ [—2, —1]. 

The choice of these parameters corresponds to have most of 
the simulated bursts with low Fo factors and with small 6 jet values. 
One could think that such distributions are already excluded by the 
observed distributions of ^jct and Fo (which are log-normal) dis- 
cussed in §1. However, those are the observed distributions of Fo 
and ^jet and they are subject to several biases (see §1). The intrinsic 
distributions might well be completely different and this motivates 
to start with this simplest assumption, i.e. that both Fo and Ojet have 
power law distributions. 

Under the hypothesis that both Fq and 9 jet have power law dis- 
tributions (with free parameters a and c varied in the above ranges 
with a step 0.2 in both parameters), only in 1% of 1000 repeated 
simulations we can find an agreement with all our observational 
constraints. In order to show the inconsistency of the simulations 
with the observational constraints we present in Fig. 2 the results 
of the simulations assuming that ^jet and Fo have power law dis- 
tributions with a = c = —1. This case, shown as an example, 
corresponds to a uniform distribution of LogFo and LogSjet. 

The rest frame Ep — Eiso plane (top left panel in Fig. 2) is filled 
uniformly with simulated bursts (yellow dots) distributed between 
the Ep — Ej limit and with a minimum 0jet=l° (the oblique right 
limit to the distribution of yellow dots). The simulated GRBs point- 
ing to us (PO) have preferentially large 9jet values (blue dots in the 
top left panel of Fig. 2). The simulated bursts of the Swift compar- 
ison sample (here represented by the smoothed density contours'^ 
- red solid lines in Fig. 2) are inconsistent with the real GRBs of 
the Swift complete sample (open squares). The red contours extend 
at high Ep values where there is a deficit of Swift bursts and they 
also over predict the number of bursts on the right hand side of the 
distribution of the real Swift bursts (i.e towards large values of E'iso 
for intermediate/high values of Ep). 

Also in the observer frame -Bpclk ~ P plane (top right panel in 
Fig. 2) the simulated Swift comparison sample (solid contours) are 
inconsistent with the real Swift bursts of the complete sample (open 
squares). Simulated bursts of the Swift comparison sample tend to 
concentrate towards the upper part of the -Epclk ^ F plane. 

The bottom panels of Fig. 2 show the cumulative rate distribu- 
tion of the fluence for the GBM and BATSE sample (right and left 
panels of Fig. 2) compared with the predictions of the simulations 
(dashed regions in the bottom panels of Fig. 2). The rate of GBM 
and BATSE bursts predicted by the simulation which assumes a 
power law distribution for both Ojet and Fq (with index -1) is a fac- 
tor ~2 larger than the rate of GBM bursts. Also the distributions 



These are obtained by staking 1000 simulations and smoothing the ob- 
tained distribution in the Ep — Eiso plane. 
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Figure 2. Result of the simulation assuming a power law distributions of ^j^t and Fq. Top left: rest frame Ep — -Eiso plane (regions and labels as in Fig 1). 
The simulation assumes a = c = —1.0 (see §4.1) which correspond to uniform distributions in the logarithm of 6jet and Fq. The yellow dots represent all the 
simulated bursts, the blue dots show the PO bursts, i.e. those pointing to us. The red solid lines are the smoothed density contour (1 and 2 a confidence levels) 
of the simulated Swift comparison sample. They should reproduce the distribution of the real Swift GRBs of the complete sample (open black squares). Top 
right: observer frame £'pjjjj. — F plane where the simulated Swift bursts (solid contours) are compared with the real GRBs of the Swift complete sample (open 
square symbols). Bottom left: cumulative rate distribution of GBM real bursts (solid black line) and its Icr uncertainty (grey solid filled region) compared with 
the cumulative rate distribution of simulated GBM bursts (solid line and dashed cyan region corresponding to its Icr uncertainty). Bottom right: cumulative 
rate distribution of BATSE real bursts (solid black line) and its \a uncertainty (pink filled solid region) compared with the prediction of the simulation (dashed 
region). 



of the peak flux of the simulated BATSE and GBM samples are 
inconsistent with the real samples. 



If we assume steeper power law distributions of Sjct and To 
[e.g. (a, c) = (—2, —2)], the excess of bursts with large peak en- 
ergy (both in the rest frame and in the observer frame of Fig. 2, top 
left and right panels respectively) is reduced but the rate of simu- 
lated GBM and BATSE bursts increases becoming more inconsis- 
tent with the real rates of GRBs detected by these two instruments 
(bottom panels of Fig. 2). This result shows that all the constraints 
that we have adopted (§3) are relevant: the GBM and BATSE com- 
parison sample map the low end of the peak flux/fluence distribu- 
tion while the Swift complete sample maps the bright burst tail of 
such distributions. The bursts of the Swift complete sample, having 
their z measured, map the distribution of GRBs in the rest frame 
-Ep - -Eiso plane. 



4.2 Peaked distributions of To and ^jct 

Since we could not find agreement between the simulations which 
assume power law distributions of 6^ct and Fq and our observa- 
tional constraints, we now consider the case of peaked distributions 

of 6'jct and Fq. 

The simplest assumption is that ^jct and/or Fo are distributed 
as broken power laws. We first assumed that only ^jct or Fo have 
a broken power law distribution, while the other parameter is dis- 
tributed as a single power law. In this case we cannot find a per- 
centage of repeated simulations larger than 2% in agreement with 
our observational constraints. 

We then considered the case of a broken power law distribu- 
tion for both Sjct and Fq: 



dN ^ r 6'f,t if 6ijet ^ 
dOjct \ ^j'et if Sjot > 



© 0000 RAS, MNRAS 000, 000-000 



8 G. Ghirlanda et al. 




48 



50 



52 54 
LogE [erg] 



56 



-5 -4 
LogF [erg/cm^] 



-3 



50 : 



50 : 




Figure 3. Result of the simulation assuming two broken power law distributions of Sjct and Fq (see text for the assumed values of the distribution parameters). 
Same plots and symbols as in Fig. 2. 



div ^ r rg if To ^ r. 
dFo 1 if To > r. 

For the distribution of ^jct we consider the following parameter 
ranges: a £ [0.5, 2.0], b £ [-2.0, -5.0] and 9, G [3°, 12°]. For 
To: c e [0.5, 2.0], d G [-2.0, -5.0] and T, G [50, 120]. The free 
parameters are varied with step 0.1 for a and b and 0.5° for 6'* for 
the broken power law distribution of 9jet and with step 0.1 for c 
and d and 10 for F, for the broken power law distribution of Fq. 

We find that at most ~20% of the 1000 repeated simula- 
tions reproduce our observational constraints when {a,b,9f) = 
(0.5,-3.0,4.5°) and (c,d,F*) = (1.8,-3.5,70) with step 0.1 
for c and d and 10 for F« for the broken power law distribution 
of Fq. a lower percentage of agreement is obtained for any other 
choice of the free parameters. 

We show in Fig. 3 the results of the simulations with the above 
parameter values for the distributions of 9 jet and Fq. We note that 
a better agreement is now found between the rate of the GBM and 
BATSE bursts (bottom panels of Fig. 3) while the distribution of 
simulated bursts of the Swift comparison sample (solid contours) 
are inconsistent with the Swift bursts of the complete sample both 
in the rest frame Ep — Eiso plane (top left panel of Fig. 3) and in 
the observer frame plane -Ep^ak ~ F (top right panel of Fig. 3). 

The assumption of a characteristic value of Fo corresponds to 
concentrate GRBs around a typical value of Ep (see Eq. 2). In this 
case the narrower 9 jet distribution reduces the number of simulated 



bursts with large values of 9 jet, thus clustering the simulated GRBs 
of the PO class around the Ep — E^ limit of Fig. 3 (top left panel) 
that was found in the case of single power laws (§4.1). 

A broken power law is a simple approximation of a peaked 
distribution. The real distribution of Fo and 9 jet could have a differ- 
ent shape. We then considered the case of log-normal distributions 
for both Fo and 9 jet, with central values of the 9 jet distribution be- 
tween 3° and 12° (step 0.5°) and width between 0.3 and 0.8 (step 
0.05) and central values of Fq between 50 and 120 (step 5) and 
width between 0.2 and 0.8 (step 0.05). We find that, if 9jet has a 
log-normal distribution with a median value of 4.5° (with a disper- 
sion of 0.5) and Fo is distributed as a log-normal with median 85 
(with a dispersion of 0.45), the 40% of the 1000 repeated simula- 
tions is in agreement with all our observational constraints. 

The latter assumption, that seems to improve the consistency 
between the simulated GRB population and the observational con- 
straints, suggests that Fo and ^jct have log-normal distributions. 
However, the fact that no more than 40% of the repeated simula- 
tions can reproduce all our observational constraints, is suggesting 
that some ingredient is still missing. This is the subject of the next 
section where we study for the first time through our numerical sim- 
ulations, the possibility that there is a relation between the average 
values of 9 jet and Fq. 
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Figure 4. Simulations assuming log-normal distributions of Bj^t and To and a relation between them (see text). Symbols and labels as in Fig. 2. 



4.3 The relation between 6jct and To 

By assuming a 6jct distribution with a characteristic value as in 
§4.2, the simulated bursts in the Ep — Eiso plane cluster around a 
correlation which is linear in this plane (i.e. parallel to the Ep — 
E-y limit), while the Ep — Eiso correlation defined by the Swift 
complete sample (and similarly by the larger, incomplete, sample 
of bursts with measured redshift - see N12) has a flatter slope, i.e. 
EpOC i?i°o . In other words, for an infinitely narrow distribution of 
djct, the simulated bursts (yellow dots in Fig. 3 top left panel) would 
produce a linear Ep — Eiso correlation which is inconsistent with 
the observed Ep — Eiso correlation. This suggests that, besides the 
fact that 6jct and To should have characteristic values (i.e. peaked 
distributions), they should also be correlated. 

Indeed, G12 find that the comoving frame properties of GRBs 
(and in particular the fact that BpOcFo and iSisoOcro^) can be com- 
bined to explain both the Epcx and the Ep(xE-y correlation 
if the ansatz ^jet^ro=const is valid. Several recent numerical sim- 
ulations of jet acceleration in GRBs suggest that a link between To 
and 9jet should exist, although the form of this relation depends 
on several assumptions of these simulations. In this section we ex- 
plore, for the first time, if a relation djet'^ro=K can account for 
the observational constraints described in §3 and in this case we 
constrain its free parameters (m and K). We start froin the result 
of the previous section, which showed that the best result (i.e. 40% 
of the repeated simulations are in agreement with the observations) 
is obtained assuming two log-normal distributions for To and Sjot- 



We simulate bursts with LogFo distributed as a Gaussians 
with a characteristic central value LogF, and a dispersion (TLogro ■ 
Similarly we assume a Gaussian distribution for Log6ljct centered 
at Log6'*jct and with a dispersion (TLogfl et- ^^^^ assume that 
there is a relation between 6'jot and Fo of the form Log6'*jct ~ 
— 1/mLogF, + q. In this way the distribution of Log6'jot is cen- 
tered on a value which is given by the assumed relation between 
6'jct and Fq. 

We explored the parameter space (defined by 5 free parame- 
ters) and found that 80% of our simulations are consistent with our 
constraints if LogF, = 1.95 with a dispersion of (TLogro ~ 0-65 
dex, m — 2.5, q = 1.45 and aLogfljet ~ 0-3 'l^''- 

We show in Fig. 4 the results of this simulation which as- 
sumes log-normal distributions of Fo and 6'jct and a relation be- 
tween these two parameters. In the Ep — Eiso plane (top left panel 
in Fig. 4) and in the -Epolk ~ F plane (top right in Fig. 4) we find 
a good agreement between the simulated Swift comparison sample 
(solid contours) and the real Swift complete sample (open squares). 
Now the predicted rate of GBM and BATSE bursts is fully con- 
sistent with the real ones (bottom left and right panels in Fig. 4 
respectively). 

We stress that, given the assumptions of our simulation (e.g. 
the spectrum, duration and unique values of the comoving frame 
energetics of all GRBs) we do not expect to find 100% of the simu- 
lations reproducing our constraints. However, we can use our code 
to derive interesting properties of the population of GRBs. Indeed, 
in our simulations we generate a population of GRBs pointing in 
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Figure 5. Distribution of Sjct of GRBs. The 9jot distribution of the total 
sample of simulated GRBs is shown by the filled circles. The solid grey 
line shows the fit with a lognormal function (Eq.7). The subsample of GRBs 
pointing towards the Earth (PO) is shown by the open (blue) squares and its 
fit with a lognormal by the cyan line. The sample of PO GRBs and peak 
flux P Js2.6 cm"'^ s~^ (i.e. the Swift comparison sample) is shown by the 
open (red) circles and its lognormal fit by the orange line. The dashed (grey) 
line shows the lognormal fit of the distribution of all the bursts (solid grey 
line) multiplied by 1 — cos 0jet. The green triangles show the distribution 
of the 27 GRBs with measured jet opening angle collected in Ghirlanda et 
al. (2004, 2007). 

every direction. Only those pointing towards the Earth (PO) are 
then compared with existing samples of GRBs (like those described 
in §3). This is also the population of bursts that will be explored by 
future GRB detectors with better sensitivity than the present ones. 
We can derive the properties of the whole GRB population (i.e. all 
the bursts pointing in whatever direction), like the jet opening angle 
distribution, the bulk Lorentz factor distribution and the true GRB 
rate. 



The GRBs that point to the Earth (PO - shown by the open 
blue squares in Fig. 5) have a Sjct distribution peaking at consider- 
ably larger values (40° - see Tab.l) than the entire GRB population. 
This can be easily interpreted: consider the distribution of the entire 
population of GRBs (black dots in Fig. 5) which contains all bursts 
pointing in every direction. The probability that a burst with a cer- 
tain 6 jet is pointing to us is proportional to ( 1 — cos ^jct ) • Therefore 
the distribution of Sjct for PO bursts is obtained from the total dis- 
tribution by multiplying by (1 — cos Sjct)- This reduces the number 
of bursts per unit 6jc% and also shifts the peak of the PO distribu- 
tion towards an average larger value. This is shown in Fig. 5 by the 
dashed (grey) line which is obtained by multiplying the fit of the 
distribution of 6jct of the entire GRB population (solid gray line in 
Fig. 5) by (1 — cos 6 jet) and it fits the distribution of the PO bursts 
(open squares in Fig. 5). 

Among the simulated bursts that are pointing towards the 
Earth we considered the bright bursts (i.e. selected with the same 
peak flux threshold of the Swift complete sample). These bursts 
tend to have small jet opening angles and this accounts for their 
9jct distribution peaking at ~5° in Fig. 5 (open red circles). 

Although apparently there is a similarity between the Oj^t dis- 
tribution of all bursts (i.e. pointing in every direction) and the Ojct 
distribution of the PO bright bursts, they differ by a factor 2 (1.8) 
in their peak values (and dispersions) which are reported in Tab. 1 . 

The three distributions shown in Fig. 5 allow us to make some 
further considerations. If we could measure ^jet for all bursts point- 
ing towards the Earth (PO in Tab. 1), we would obtain the open 
(blue) square distribution of Fig. 5 with a mean ~ 40°. However, 
the real Ojct distribution of the population of GRBs (i.e. all the sim- 
ulated bursts - black filled circles distribution in Fig. 5) has a mean 
of ~8.7° and it is more consistent with the distribution of the sim- 
ulated PO bursts with large peak fluxes (the Swift comparison sam- 
ple). This suggests that the bursts distributed in the high part of the 
Ep — Eiso correlation, where are the bursts of the complete Swift 
sample (filled black dots in Fig. 4 top left panel), properly sample 
the peak of the 6jct distribution of the entire GRB population. 



4.4 6jct distribution of GRBs 

From the best simulation described in §4.3 we can derive the dis- 
tribution of the jet opening angle of GRBs. In Fig. 5 we show the 
distribution of 6 jet for all the simulated bursts (black points) and for 
the PO bursts (open cyan squares). The population of GRBs point- 
ing towards the Earth and with a peak flux P Jj2.6 cm^^ s^^ in the 
15-150 keV range (i.e. the Swift comparison sample) is shown by 
the open (red) circles. All the distributions of ^jet can be modeled 
with a log normal function: 



N(x) 



A 



xa\'2Tr 



exp 



{Inx — 



2cr2 



(7) 



where the free parameters are (/i, a) and the normalization A. The 
best fit parameters fi and a are reported in Tab. 1 . The peak of the 
log-normal distribution, i.e. its mode, is exp{fj, — cr^), the mean is 
exp(/i + 0-^/2) and the median is exp{ff). Since the asymmetry of 
the log-normal distributions can be considerably large, we report 
in Tab. 1 all these moments. 

The 6jct of GRBs of the Swift comparison sample (red open 
circles in Fig. 5) have a mean of Sjct~ 4.7°. This distribution is 
consistent with the 9 jet estimated from the break of the optical light 
curves (Ghirlanda et al. 2004, 2007), shown by the open (green) 
triangles in Fig. 5. 



4.5 GRBs with no jet break 

It has been shown in §3 that if a burst has a To such that sin 6jct ^ 
l/Fo, its Eiso is determined by To (Eq. 5) and not by 6 jet- This 
value is lower than that computed by 6jet (Eq. 4). In these bursts, 
therefore, we should not observe a jet break in their light curve 
since the emitted radiation is initially collimated within an angle 
arcsinl/Fo larger than ^jet. Since F decreases during the after- 
glow phase due to the deceleration of the fireball by the interstellar 
medium, in these bursts the jet break, corresponding to the transi- 
tion l/F ~6jct, will never happen. 

The above argument contributes to explain the fact that bursts 
might not show an evident jet break in their afterglow light curve 
if l/Fo^ sinfljot. However, in these bursts we expect that the after- 
glow light curve is declining with a typical post-break decay index 
~ —p (where p is the shock-accelerated electron energy distribu- 
tion index - e.g. Panaitescu & Kumar 2001). Other possible expla- 
nations for the lack of fbreak measurements have been proposed. 
Numerical simulations (e.g. Van Eerten et al. 2010), for instance, 
suggest that the jet break transition can be very smooth (almost 
difficult to be distinguished from a single power law decay with 
available data sets) due to a combination of the jet dynamics be- 
fore and after the jet break time (and additional complications can 
be induced by the viewing angle effects when the observer is not 
on-axis). Although a detailed discussion of the missing jet breaks 
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Table 1. Parameter values (fj, and a) obtained by fitting a log-normal function (Eq. 7) to the distributions of 8 jet (Fig- 5, 6) and Fq (Fig. 7) for all the simulated 
bursts (ALL), for those pointing to us (PO) and for those pointing to us and with a peak flux larger than Piim, i.e. the Swift comparison sample (PO Swift). (*) 
fit of the distributions of GRBs pointing towards the Earth that should not have a jet break (see §5.2). For each distribution are reported the three moments: the 
mode, the mean and the median. 
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Figure 6. Distribution of 9jct of GRBs: the PO simulated GRBs (open blue 
squares) and the GRBs of the PO class that have sin Bj^t ^ 1 /To (open 
green triangles) are shown. The latter are those that should not show any jet 
break time in their afterglow light curve. For the PO bursts with P ^2.6 
cm~^ s~^ (open red circles) we show the subsample of bursts that have 
sinSjct l/Fo (open orange stars). 



in GRBs is out of the scope of this paper, we notice that bursts with 
sinSjet ^ l/Fo can partly account for the explanation of the lack 
of measured jet breaks. This is the first time that such an argument 
is presented and surely deserves further studies. 

Fig. 6 shows the distribution 6jct of PO bursts (open blue 
squares) and the subsample of bursts with no jet break (open green 
triangles). These amount to ~6% of PO bursts. The mean of their 
log-normal distribution is 6jct~20°. One testable observational 
prediction of our simulations is that GRBs with no jet breaks should 
be preferentially soft (E°^^ of few tens of keV) The open red circles 
in Fig. 6 correspond to PO bursts of the Swift comparison sample 
while the open orange star symbols correspond to bursts with no 
jet break. These have a mean jet opening angle ~3.3°. We find that 
~2% of the Swift bright bursts should not have jet break in their af- 
terglow light curves. They could correspond to those events which 
do not show any evidence of a jet break in their optical light curve 
(e.g. Mundell et al. 2006; Grupe et al. 2007) although other obser- 
vational selection effects very likely contribute to the paucity of the 
jet break measurements. The fit of the distributions shown in Fig. 6 
with log-normal functions are reported in Tab. 1 . 



Figure 7. Distiibution of Fq of GRBs. Symbols as in Fig.5. The 30 GRBs 
with Fo estimated from the peak of their afterglow light curves (G12) are 
shown with the green open triangles. 



4.6 To distribution of GRBs 

From our simulation we can derive the distribution of To (Fig. 7). 
The total population of simulated bursts (filled circles in Fig. 7) 
has a log normal distribution with a mean To =274. Those pointing 
towards the Earth (open blue squares in Fig. 7) have a smaller mean, 
ro=49. The PO bursts with peak flux larger than 2.6 cm~^ s~^, 
i.e. those of the Swift comparison sample, have a typical ro=355. 
Although the distribution of Fq factors for those bursts with a peak 
in their afterglow light curves (G12) is still made of few events, it 
agrees (open green triangles in Fig. 7) with that predicted by our 
simulations (for the sample of PO bursts of the Swift comparison 
sample - open red circles in Fig. 7). 

Also in the case of Fo we note that if we were able to mea- 
sure Fo for all the bursts that point towards the Earth, we would 
obtain a slightly smaller peak value of Fq with respect to that of the 
distribution of all the bursts (pointing in every direction. 

We note that the Fq distribution of the general population of 
GRBs peaks at considerably low values of Fq. This is a result of our 
simulations where, as explained in §4.2, we assume a peaked log- 
arithmic distribution of Fq with free peak and width. If we assume 
a distribution of Fq with a smaller fraction of bursts with low Fq- 
values, then we cannot reproduce the flux and fluence distributions 
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Figure 8. GRB rate as a function of redshift. The sample of simulated bursts 
is shown by the filled black circles and the GRB formation rate assumed in 
our simulation (see §2) is shown by the solid grey Une. This curve is normal- 
ized to the histogram. The GRB formation rate without density evolution (as 
derived by Li et al. 2010 - see Eq. 1) is shown by the dashed grey line. The 
rate of bursts pointing to us (PO) is shown by the open (blue) squares and 
that of the PO bursts with P ^2.6 cm~^ s~^ (i.e. the Swift comparison 
sample - PO Swift) by the open (red) circles. For comparison is also shown 
the rate of Swift GRBs of the complete sample (rescaled to match the rate of 
PO Swift rate). The solid lines reported in the top of the plot are the cosmic 
rates of SNIb/c computed by Grieco et al. (2012) for different assumptions 
on the cosmic star formation rate. 

and the detection rates of GRB detection of the GBM and BATSE 
instruments. Therefore, our simulations predict that a considerable 
fraction of GRBs should have To as low as a few tens. These bursts 
might well be detected by current instruments. While the detailed 
study of their prompt and afterglow properties is out of the scope of 
the present paper, we note that their prompt emission should hardly 
differ from that of bursts with larger Fq values (except for the ob- 
vious fact that their prompt and Liso is lower). In fact, if To 
is low the fireball deceleration timescale (e.g. Eq. 14 in Ghirlanda 
et al. 2011) is tpcak ~ '^^llo50^o\^^ hours which is much larger 
than the prompt emission timescale. So, while the prompt emission 
of low-Fo burst should not be influenced by the afterglow contri- 
bution, their late time afterglow onset could be a distinctive feature 
(typical afterglow onset timescales are of the order of few hundreds 
second - Ghirlanda et al. 201 1). 



4.7 The GRB rate 

Another consequence of our simulations is the rate of GRBs. This 
is shown as a function of redshift for the entire population of simu- 
lated bursts (filled circles in Fig. 8), in units of bursts Gpc~^ yr~^. 
The GRB redshift distribution (Eq. 1) assumed in our simulations 
is shown by the solid grey line in Fig. 8 and the observed star for- 
mation rate (Li 2008) is shown by the dashed (grey) line rescaled 
by an arbitrary factor to match the rate of GRBs at z = 0. We 
also show the rate of PO bursts (open blue squares) and that of the 
PO bursts of the Swift comparison sample (open red circles). Fig. 
8 also shows the recent estimate of the rate of SNIb/c computed by 
Grieco et al. (2012). The different curves for SNIb/c correspond to 
different assumption of the cosmic star formation rate (CSFR) in 
that paper. As a result of our simulation, the local rate of GRBs is 
~0.3% that of SNIb/c. 



5 SUMMARY AND DISCUSSION 

We have studied two fundamental parameters of GRBs: the jet 
opening angle Sjot and the bulk Lorentz factor Fo. The first ques- 
tion that we aimed to answer was whether Ojct and Fo have prefer- 
ential vahies. The direct measure of 9 jet through the jet break times 
observed in the optical light curves (Frail et al. 2001; Ghirlanda et 
al. 2004, 2007) shows that ejot~5°. The measure of Fo from the 
peak of the afterglow light curve for ~30 GRBs (G12) also shows 
a characteristic value* of Fo~60. However, the limited number of 
events with a direct estimate of Sjct and Fo and the possible se- 
lection effects, related to the difficulties of measuring these two 
parameters (see §1), prevent us to assume them as representative 
of the GRB population. In particular we want to test the consis- 
tency of different possible distributions of 6 jet and Fo with a set of 
available observational constraints (§2). Moreover, we aim at con- 
straining the free parameters of the distributions of Oj^t and Fo and 
derive if and how these two parameters are correlated. 

In this paper we used a population synthesis code to simulate 
GRBs with different assigned distributions of Ojct and of Fo each 
one with a set of free parameters that we left free to vary within 
certain ranges. Obviously, we did not assume the observed distri- 
butions of ^jct and Fo as constraints to avoid circularity. 

We assume that GRBs have a unique comoving frame peak 
energy E'^ and collimation-corrected energy E'^ (the large black 
dot in Fig. 1) which are transformed into their corresponding rest 
frame E^ and E^ respectively. The assigned 9jct and Fq allow us 
to derive the isotropic equivalent energy of the simulated bursts 
according to the relative value of ^jct and Fq. Eiso^E-y /O'^^^ if 
sin6'jot > I/Fq, while Eiso'-^E-^V'Q in the opposite case. This in- 
troduces a "natural bias" in the distribution of Eiso'- those bursts 
with a small enough Fo will have an isotropic energy which is 
smaller than that one would calculate using the value of djet- In the 

1/3 

Ep — Eiso plane of Fig. 1 this corresponds to a limit EpCC E-J^ . 
Bursts with 1/Fo>sin6'jot will lie along this limiting line and there 
should be no GRBs on the right of this line [i.e. in region (III) in 
Fig. 1]. 

This is the first time that the limit mentioned above is consid- 
ered within the framework of studying the distributions of GRBs 
in e.g. the Ep — Eiso plane. Indeed, this limiting line can account 
for the absence, in the observed GRB sample with measured z and 
well constrained peak energy (i.e. the bursts used to construct the 
Ep — Eiso correlation), of bursts with intermediate/low peak energy 
Ep and very large iSiso- 

The assumed distributions of Fo and ^jct determine the dis- 
tribution of simulated bursts in the Ep — Eiso plane of Fig. 1 . We 
considered two types of distributions for ^jct and Fq: (A) a power 
law distribution, i.e. 6'jot and Fq do not assume any preferential 
value or (B) both ^jct and Fo have peaked distributions (either bro- 
ken power law or a log-normal distributions). 

In order to test these two hypothesis we compared the results 
of our simulations with three GRB samples: the complete Swift 
sample of GRBs detected by BAT with measured redshifts (S12), 
the sample of bursts detected by the GBM in the last 2 years (Gold- 
stein et al. 2012) and the 4th BATSE catalog of GRBs (Meegan et 
al. 1997). The simulations should reproduce several proprieties of 
these samples. 

While most of the bright bursts of the Swift complete sample 
of S12 have measured z and provide an observational constrain in 

* This average value is obtained assuming that the circumburst medium 
has a wind density profile (see G12). 
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the rest frame — Eiso and observer frame -EpoJji — F plane (left 
and right panels of Fig. 2,3,4, respectively), the number count dis- 
tribution and rate of the BATSE and GBM populations of bursts 
(mostly without measured z) are used as additional constraints 
since they map the faint end of the number count distribution of 
GRBs. 

Our main result is that we cannot reproduce all our observa- 
tional constraints if the To and Ojet distributions are power laws. In 
this case the rate of GBM and BATSE bursts predicted by our sim- 
ulations is a factor ~2 larger than the real one and the distribution 
of the Swift simulated bursts in the Ep ~ EiBo and E°^ly. — F plane 
is inconsistent with the real complete sample of Swift bursts. 

Instead, if ^jet and To have broken power law distributions 
(with peak values 6'jot~ 4.5° and r()~ 70) or log-normal distribu- 
tions (with peak values ^jct~ 4.5° and Fo--^ 85) a better agreement 
between the simulations and the observational constraints is found. 
However, the broken power law or log-normal case produce a lin- 
ear Ep — Eiso correlation due to the assumption that the simulated 
bursts have a ^jct distribution with a unique peak value (see §4). 
This motivated us to consider the possibility that there is a relation 
between the peak values of the distributions of ^jet and Fo. G12 
foimd that among GRBs with a Fo estimate, three new correlations 
are found: SisoOcFo^, LisoOcFo^ and SpOcFo. The combination of 
these correlations with the assumptions that (?j?c5or()=const allows to 
derive the three main empirical correlations of GRBs: the Ep—Eiso 
correlation, the Ep — Liso correlation and the Ep — E^ correlation. 

We therefore assumed that both Ojet and Fo have log-normal 
distributions and that a relation of the type 0j2;Fo=const exists be- 
tween the peak values of their respective log-normal distributions. 
We found good consistency between our simulations and the obser- 
vational constraints (Fig. 4) in the case of a log-normal distribution 
of F() with central value 90 and logarithmic dispersion of 0.65. The 
distribution of ^jct (also a log-normal) is in this case determined by 
the relation ^j^tFo=const which we find should have m = 2.5. This 
value is what one obtains by combining the above scaling relations 
(between Fo and Eiso, Ep) with the Ep — Eiso correlation of the 
Swift complete sample which is Epcx Efso ■ The existence of a rela- 
tion 6'j"\r()=const (with rn ^ 1 and const= 10-40) is also predicted 
from recent models of magnetically accelerated jets in GRBs (e.g. 
Tchekhovskoy et al. 201 1). 

We found that the Fo distribution that best reproduces all our 
observational constraints extends to low Fo values. If wc cut the 
Fo distribution so to exclude such low values of Fo we cannot re- 
produce the observed fiux and fluence distributions and detection 
rates of BATSE and GBM. Therefore, we find that low-Fo bursts 
should exists in populations of GRBs detected by most sensitive 
detectors. Although a detailed study of the prompt and afterglow 
properties of these events is out of the scopes of this paper, we can 
draw some remarks. Apart from their relatively low Ep and Liso 
(which are correlated with Fo as found by Gil), the low-Fo bursts 
should have a late time afterglow onset (i.e. a few hours for typical 
parameters, see §4.6). Therefore, their prompt emission should not 
be contaminated by the afterglow while their late time afterglow 
onset could be one of their distinctive features. 

An immediate consequence of our results is that the large scat- 
ter of the Ep — Eiso correlation can be interpreted as due to the jet 
opening angle distribution of GRBs. The found inverse relation be- 
tween Sjet and Fo implies that bursts with the largest bulk Lorentz 
factors should have a smaller average Ojet- On the other hand, bursts 
with relatively low average Fo factors should also have, on average, 
large 6'jet. 

Our results depend on the assumption that all bursts have the 



same E'p=l.5 keV and £'^,=1.5x 10"'* erg. Although there could 
be a dispersion of these values, our results still hold if the width 
of this dispersion is not larger than the dispersion of the observed 
quantities. We note that larger values of Ep and E!y would move the 
EpxEiso^''^ (Eq. 6) towards the upper part of the plane of Fig. 1. 
As a consequence some of the real GRBs of the Swift complete 
sample, would be cut out of the plane because they would lie in the 
forbidden region (in) of this plane. On the other side we could as- 
sume lower values of Ep and £"'.. Since we do not know their real 
dispersion, we tried to assume £'p=0. 15 keV and E!y=l.5xl0'^'^ 
(i.e. a factor 10 lower than the values assumed in the simulation). 
Under this different assumption, for the case of log-normal dis- 
tributions of both Ojet and Fo and of an intrinsic relation between 
these two parameters, we find that the Ojet distribution is consistent 
with that found with the fiducial values of Ep and Eiy, but with 
a different distribution of Fo. Indeed, in this case we find a mean 
value Fo~ a factor 3 larger than that of the present simulation. Al- 
though it is not possible at the present stage to constrain the dis- 
tribution of Ep and E'^ , these results suggest that their dispersion 
should be lower than a factor of ~10. 

Our best simulations allow us to derive the properties of three 
populations of GRBs: those that are pointing to us and that have 
a peak flux bright enough to enter in the Swift bright sample (i.e. 
with the same peak flux threshold adopted for the Swift complete 
sample of S12), those that are pointing to us and, finally the full 
population of simulated GRBs, oriented randomly in the Universe 
(i.e. pointing to us and not). The latter is the GRB population that 
we cannot study on the base of the bursts that we detect. The main 
advantage of our population synthesis code is that we can infer the 
properties (e.g. the Fo and Ojet distribution and the true GRB rate 
in this work) of this population of bursts, which is unaccessible 
through the observations. 

One immediate consequence of our simulation is the true 

-Ep — Ei^o correlation. If we consider the PO bursts and if we 
were in principle able to detect them all, we should find a differ- 
ent Ep ~ Eiso correlation than the one presently reported in the 
literatiu'e. Indeed, the fit of the PO bursts in the Ep — Eiso plane of 
Fig. 4 yields a correlation with slope 0.5 and normalization -27.6 
while the entire GRB population, the total simulated bursts, have 
a Ep — Eiso correlation with slope 0.44 and normalization -20.7. 
This is due to the fact that PO bursts tend to populate the lower re- 

1/3 

gion of the Ep — Eiso plane (Fig.4) where the E^J^ limit cuts their 
distribution in the Ep — Eiso plane. Therefore, if we could measure 
-Eiso and Ep for all the bursts that point to us, we should determine 
a flatter Ep — Eiso correlation than that observed so far in the high 
part of the plane with bright bursts. 

Our simulation predicts that the bright bursts detected by Swift 
should have a mean opening angle of ^jet~4.7°. This value is only 
a factor 2 smaller than the mean of the entire GRB population that 
we have simulated (which has Ojet^ 8.7°). However, from Fig. 
5 (open blue squares) one can see that if we were able to detect 
fainter GRBs and to measure their jet opening angle, we would 
obtain a mean of 40° . Intriguingly we note that the present distri- 
bution of ^jct measured from the optical afterglow break times in a 
few bursts is representative of the Ojet distribution of the entire pop- 
ulation of bursts. This is because the bursts that we have detected 
so far populate the high region of the Ep — Eiso plane where the 
Ojet distribution can be almost unbiasedly sampled. In fact, only the 
bursts at lower values of Ep and Eiso are affected by the "natural 
bias" of l/Fo>sin^jet. The low Ep- low Eiso region is where PO 
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bursts concentrate (they have large ^jet or small To, enhancing the 
probability to point at us). 

Our simulation predicts that there are bursts with no jet break, 
the ones with l/ro>sin^jet. Their afterglows will never have a jet 
break since the condition l/To^sin^jet is never met but their af- 
terglow light curve should have a characteristic post-jet break in- 
termediate/steep decay slope. These should be ~6% of the bursts 
pointing to us and ~2% of the bursts detected by Swift with 

P>flim. 

According to our best simulation, the mean To of all bursts is 
(Fo) = 274. The To distribution is highly asymmetric and there is 
a considerable difference between its mode (i.e. the peak) the mean 
and the median. The simulated bursts pointing to us corresponding 
to the Swift complete sample have (Fo) = 355. These two values 
are broadly consistent, as explained above, since these bursts pop- 
ulate the upper part of the — E\so plane where the distribution 
of GRBs is almost free from the "natural bias". Remarkably, if we 
were able to measure Fo for all the bursts pointing to us, we would 
find a very low value of the mean of (Fo) = 50. Finally, we have 
found that the distribution of Fo that we predict for the Swift bright 
sample is consistent with the distribution of Fo of the GRBs studied 
in G12. 

We can derive from our simulations the true rate of GRBs. 
Previous studies of the GRB rate assumed a unique value of 6'jot, 
typically 0.2 rad or the observed distribution of 6'jet(e.g. Guetta et 
al. 2005; Grieco et al. 2012). Our simulations (§4.5) show that the 
peak of the intrinsic/global distribution of ^jet is a factor 2 larger 
than the real intrinsic distribution and has a much wider dispersion 
(Tab. 1). Differently from existing GRB rate estimates based on the 
correction of the isotropic GRB rate for an average beaming fac- 
tor (e.g. Guetta et al. 2005; Grieco et al. 2012) in our simulation 
the total number of simulated bursts is adjusted in order to repro- 
duce the rate of detections of GBM and BATSE. Therefore, we 
have the rate of GRBs as a function of redshift independently from 
the value of 6 jet of each single burst. If we compare this rate with 
that of SNIb/c (from Gireco et al. 2012) we find that the local rate 
of GRBs is ~0.3%. Moreover, if we consider the 7% fraction of 
SNIb/c which produce Hypernovae events (Guetta & Delia Valle 
2007) we find that the rate about 4.3% of local Hypernovae should 
produce a GRB. 
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